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In fl| Willis and Kerswell study in direct numerical 
simulations the statistics of turbulent lifetimes in pipe 
flow. Their simulations of a very long pipe are of consid- 
erable interest in relation to the chaotic saddle picture of 
the transition to turbulence 2]. They suggest that their 
data for six different Reynolds numbers support a diver- 
gence of the lifetime near a Reynolds number of about 
1870. However, their conclusion is not compelling: a re- 
analysis of their data shows that it is compatible with an 
exponentially increasing lifetime, as observed in [5]. 

The results in Jj are based on the turbulent lifetimes 
of flows which are prepared from snapshots of a turbu- 
lent simulation at Reynolds number 1900 which is then 
integrated at a lower Reynolds number until they relam- 
inarize or until they reach a maximal integration time 
of 1000 U/D. From 40 to 60 such initial conditions the 
probability P{T) to remain turbulent for a time T at 
least is obtained and compared to an exponential varia- 
tion, P{T) ~ exp(— Tln2/r/i), with Th{Re) the charac- 
teristic lifetime. The problematic part is that this expo- 
nential distribution is expected to occur for long times 
only and that the choice of preparation of initial con- 
ditions introduces transient elements into P{t) of unpre- 
dictable duration. This is evident from a magnification 
of their Fig. 3, shown here as Fig. 1. P{T) for the lowest 
Reynolds number Re = 1580 clearly is constant up to 
about 38 U/D, drops off slowly until 58 U/D and then 
falls of more steeply for longer times. An exponential fit 
to the tail for times T > 58 gives a lifetime of Th ~ 4.9, 
much shorter than the value that can be read off from 
their Fig. 2, Tfi,wK ~ 14.5. The lifetime statistics at 
Re = 1700, 1740, 1780, and 1820 also show several re- 
gions with different behaviour, and again only the tails 
should be fitted. The data at Re — 1860 is inconclusive, 
since it covers a small range of P{T) only. 

The re-analysis of the slopes has consequences for the 
variations of lifetimes with Reynolds number, as shown 
in Fig. 2. The reduced lifetimes at Re = 1580 give 
a higher value for l/rh- The lifetimes extracted from 
Fig. 1 lie close to the experimental data of 3], and are 
compatible with an exponential variation with Reynolds 
number. A logarithmic representation shows that the 
last data point lies somewhat lower, and that there 
are small differences in the slopes. Unfortunately, the 
simulations stop where the experiments begin, so it is 
difficult to see whether the observed differences point to 
a transition in the character of the turbulent dynamics. 
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FIG. 1; Turbulent lifetimes P{T) for pipe flow. Shown is 
a magnification of Fig. 3 from l| (symbols and Reynolds 
numbers as used there) together with straight lines indicating 
the fit to an asymptotic exponential (straight line). The 
dashed line for Re — 1580 has the slope they use in their 
analysis. 








1 

b 

0.1 
0.01 
0.001 

- 




^ 0.0001 

\ . \ . ^ 



+ 

1600 1800 2000 







1 





1600 1800 2000 

Re 

FIG. 2: Variation of the half lifetime with Re. The red 
crosses are the data of the circles the ones read off from 
the fits shown in Fig. 1 and the margenta triangles are 
the data of 0]. The dashed line is a fit to an exponential 
variation. The inset shows the same data on a logarithmic 
scale. 


We conclude that the reanalysis of the data shows that 
they are compatible with an exponential increase of the 
lifetimes, as observed in 
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